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Abstract 

Motivated in part by string theory, we consider a modification of the ACDM cosmological 
model in which the dark matter has a long-range scalar force screened by light particles. 
Scalar forces can have interesting effects on structure formation: the main example presented 
here is the expulsion of dark matter halos from low density regions, or voids, in the galaxy 
distribution. 



February 1, 2008 



*e-mail: ssgubser@Princeton.EDU 
t e-mail: pjep@Princeton.EDU 



1 Introduction 



What if there were a long-range force other than gravity acting only on the dark matter? 
The idea of a long-range force from the exchange of a massless scalar has a long history, with 
more recent attention on a scalar force in the dark sector that avoids the constraint from the 
Eotvos experiment (as reviewed in pQ). A screening length on a dark sector force generated 
dynamically by the presence of light particles with a Yukawa coupling to the scalar was 
discussed in £Q , and the picture of two dark matter species with a scalar force of attraction 
among like particles and repulsion among unlike particles was discussed in [2]. Here we 
combine the two notions, in the dark matter lagrangian 

£ = \{ d< t>) 2 + + <My#+ + t_iy*_ 

- y s <p^ s ^s - (m+ + - (m_ - . 

The constants m± and y± are all positive. The fermions ^± are the non-relativistic dark 
matter, and the quanta of the additional species ^> s will be termed screening particles. The 
screening particles are massless at = 0, and they can provide dynamical screening by 
pulling the field close to <fi = 0. An advantage of this mechanism is that the screening length 
increases with the size of the universe, so that its length in co-moving coordinates is roughly 
constant: thus it can play an important role in the formation and evolution of the galaxies 
at low redshift but not substantially affect the formation of the observed anisotropy of the 
3 K thermal background radiation (the CMB) at decoupling at redshift z ~ 1000. 

We have two purposes in this continuation of our previous work QjJ First, we con- 
tinue the exploration of possible astrophysical implications of the scalar-mediated force that 
emerges from the action (fT|). in this paper with particular attention to the properties of 
voids. We argue that the repulsive force between dark matter particles with opposite scalar 
charge can make the voids more empty than the prediction of the ACDM model, and an 
apparently better approximation to what is observed [3]. Second, we argue that the ac- 
tion fits rather neatly into the framework of supersymmetry and string theory, although 
supersymmetry breaking as it is usually understood causes some problems. In summary, we 
argue that a scalar-mediated, dynamically screened interaction in the dark sector has both 
theoretical motivation and possibly interesting observational consequences. 

In the modified CDM model in |2j, the dark matter could be described in terms of a 
scalar field coupled to two species of fermions, with an interaction lagrangian of the form 

C mt = -moe^"^!?! - m e- WM ^ 2 ^2 , (2) 

where m is close to the Planck scale and (3 is of order unity. The scalar interacts only with 
these particles, giving rise to an inverse square law proportional to the product of "scalar 
charges" Qi = drrii/difi. The total scalar charge density must vanish once <fi relaxes to its 
equilibrium point, determined by an effective potential that includes contributions from the 
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Figure 1: The solid lines depict the masses for non-relativistic dark matter species included 
in the model of |2j. The dashed line is m(0) for the additional, light species while the 
gray line shows the physical mass |m(</>)|. Including all three species leads to a model like 
(UJ, provided <j) stays close to the zero of m(4>). 

dark matter. Additional species with couplings to the scalar do not change the picture much, 
provided they are heavy. 

The variant (JTJ which we explore in this paper is based on the idea that the mass 
parameter for some species could pass through zero. For example, the mass term for \l/ s is 
— Hs^s^s, but the overall sign of this term can be reversed by performing a chiral rotation, 
ty s — > e lnl5 fy s . Thus for (^J), the mass of ^ s is \y s <j>\- The situation is summarized in figure^ 
If there are many \l/ s quanta compared to the number of ^f± quanta, then <fi is pinned near 
zero by the contribution of \l/ s to the effective potential. The total scalar charge for ty± 
quanta need not be zero because the screening particles can make up the difference — a first 
hint that effectively screens the scalar interaction. 

The organization of the remainder of the paper is as follows. In section |2] we describe some 
reasons to regard scalar forces in the dark sector as well-motivated in string theory, and we 
construct supersymmetric generalizations and extensions of the action (fT|)from ingredients 
that arise in compactifications of string theory. 1 In section |3] we briefly recapitulate the 
dynamical screening mechanism. Aspects of the growth of structure in the dark matter 
distribution, in linear perturbation theory and a spherical model, are discussed in section EJ 
with particular attention to the case where the two dark matter species have very different 
scalar charge-to-mass ratios. In section 14.31 we consider the possible application to the 
properties of voids in the distribution of normal galaxies. 

2 Supersymmetric and string theoretic models 

There are both top-down and bottom-up motivations to consider a link between string theory 
and scalar forces in the dark sector. 

lr rhe reader wishing to pass lightly over the string theory may find the summary in subsection 12 . 51 helpful . 
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A top-down motivation is that there are moduli in string theory, and there are heavy ob- 
jects whose mass depends on them, offering the opportunity to stabilize moduli not through 
vacuum effects but through the presence of some density of the heavy objects jHE]- It is hard 
to see how this could work in the visible sector because of the tight constraints on fifth forces 
and on the variations of fundamental constants such as the fine-structure constant. But in 
the dark sector, it is not only possible — because of the considerably weaker constraints on 
the Equivalence Principle for dark matter particles — but interesting for structure formation, 
as has been argued in QjJ |2| and will be further argued in section 14.81 

A bottom- up motivation is that, to have interesting effects, scalars in the dark sector must 
be very light, and the only obvious means to accomplish this is supersymmetry. (Below we 
will review the situation when supersymmetry is broken). More specifically, with the interac- 
tions in ([TJ present, it is unnatural for the potential V(4>) to be zero: absent some symmetry 
argument, the natural scale for V(<f>) is the Planck scale. So let us provisionally think of <p 
as a flat direction whose existence is protected by standard non-renormalization theorems of 
unbroken supersymmetry. 2 Actually, in a supersymmetric theory, massless scalars are com- 
plex, and there could be several of them. But the physical consequences are fairly similar to 
the simplest case of one real scalar. 

A further qualitative feature that (0) provides is degrees of freedom (the screening par- 
ticles) which become massless at a special value of the massless scalar 0. The spin and 
statistics of these particles is not important: the density in single particle phase space can 
be negligibly small, and the dark matter Compton wavelengths are very much smaller than 
length scales of interest in astronomy. But the occurrence of additional massless particles at 
special points in the moduli space is a hallmark of string compactifications, and there are 
even some dynamical reasons to think that vacua with such particles present are preferred 
[Hj. Such vacua are often referred to as having enhanced symmetry because a particularly 
common example is to have some non-abelian gauge symmetry become unbroken. In fact, 
the logic of jH] is that if the universe evolves through a configuration with extra massless 
particles, then these particles will be produced by quantum effects and will tend to prevent 
the further evolution from moving away from the point of enhanced symmetry. One of our 
desiderata, to be described in more detail in section |HJ is that the screening particles should 
be large in number compared to the massive dark matter particles, and the mechanism of 
[B] will tend to make this so. 3 

The massive dark matter particles need to have a ratio of scalar charge to mass on the 
order of 1/Mpi, though with screening, it is interesting to consider such ratios an order of 
magnitude higher. String moduli are associated with the gravitational physics of the extra 
dimensions, and the massive dark matter particles can be strings or branes stretched around 
some cycle or between some pair of other branes: then a charge-to-mass ratio on the order 

2 Another way to get a light scalar is to have it be a Goldstone boson or axion, but the Yukawa-type 
couplings in Q arc forbidden. 

3 However, as we will describe in section Ts. 31 the screening particles need to be highly energetic — more so 
than the mechanism of 6 would seem likely to predict. 
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of 1/Mpi is indeed the natural range. 

In summary, there is an interesting convergence of lines of thinking toward models of the 
type (0). In the following subsections, we will consider two super symmetric models in more 
detail, note some pitfalls, and present one further argument in favor of light scalars even after 
supersymmetry breaking. The first model, based on chiral superfields, fails to reproduce the 
physics of but for interesting reasons which will guide us to a better model, presented 
in section YIM 

2.1 A chiral model and the pitfall of D-terms 

One choice for the dynamics is to replace and \l/ s with three chiral superfields, $1, $2, and 
$3, with canonical kinetic terms and a superpotential 

W = y s ^ 2 ^ 3 . (3) 

Assuming standard kinetic terms, the lagrangian for the component fields of the $j only is 

^ = Yl (\ d( Pk\ 2 + -^fki^k) - (^01^2^3 + VsMzIpl + ^03^1^2 + h.C.) 

k V 2 1 (4) 

where we decompose a Majorana spinor as ^ = f|d) ; an d ipiip 2 means ipfip 2a - In the 
flat directions for the theory one of the 0, is non-zero and the others vanish. Suppose 
0i 7^ 0. Then the bosonic and fermionic fields of $2 and $3 acquire masses |y s 0i|, but 0i is 
protected against acquiring a potential by the supersymmetry. Thus 0i plays the role of 
in the lagrangian and the component fields of $2 and $3 play the role of 

We may reintroduce the non-relativistic dark matter species \l/± as components of two 
additional chiral superfields <3>±, again with standard kinetic terms, and with additional 
terms in the superpotential of the form 

SW = ^(m + + + - . (5) 

Instead of writing out the component lagrangian in detail, let us simply remark that 0± and 
wind up having masses m± ± y±0i, so they play the role of the fields in (Q): their 
quanta are the non-relativistic dark matter. 

In [7|, strong evidence was presented for the claim that M-theory on R 3 ' 1 times a par- 
ticular seven-manifold of G2 holonomy, asymptotic to a cone over SU(3)/U(1) 2 , has for its 
low-energy dynamics the supersymmetric model described by three chiral superfields with 
superpotential y s $i$ 2 < ^ ) 3 (see figure^). The G 2 cone is not compact: it has infinite volume. 
So gravity is decoupled and the theory has only global supersymmetry. But asymptotically 
conical G 2 manifolds describe the local geometry of compact G 2 manifolds near isolated 
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Figure 2: a) A cartoon of the G 2 manifold asymptotic to a cone over SU(3)/U(1) 2 . The 
topology of the manifold is CP 2 x R 3 . b) M-theory dynamics on the G2 manifold is equivalent 
to type IIA string theory with a three D6-branes intersecting at a point. 

pointlike singularities, and in the compact case gravity ceases to decouple and supersymme- 
try becomes local. Thus it would seem that a G2 compactification that includes a singularity 
where some SU(3)/U(1) 2 shrinks would be a good candidate for an embedding of (JTJ) in M- 
theory. The massive dark matter particles could be wrapped M2-branes that pass through 
or close to the singularity. The appeal of this picture is enhanced by work (see for exam- 
ple [HI El) indicating that G2 manifolds with singularities (more precisely, intersecting loci 
of ADE singularities) offer a way of constructing quasi-realistic theories that include the 
Standard Model. 

However, there is a significant hitch. Consider the type IIA description of the cone over 
SU{3) /U{1) 2 : three D6-branes intersecting at a point (see figure|2t>). The fields $1, $ 2 , and 
$3 are strings stretching between pairs of branes. After compactification, the gauge fields 
on the D6-branes must be expected to become dynamical. Their coupling to the fields $j, 
together with Af = I supersymmetry, implies a D-term potential of the form 

V D = g\ (|0 2 | 2 - |0 3 | 2 ) 2 + gl (|0 3 | 2 - I0i| 2 ) 2 + gl (|0i| 2 - I0 2 | 2 ) 2 , (6) 

where the gi are the gauge couplings for each D6-brane. The potential (jUJ), together with the 
F-term potential exhibited in the second line of (J3J), leaves no flat directions: for example, if 
one sets 4>2 = 4>3 = 0, then the potential in the 0i direction is (gf + 5 , l)l ( / ) i| 4 - The fact that D- 
term potentials can lead to directions where a quartic self-interaction (with coupling of order 
unity) keeps the field pinned to zero is interesting for the version of chameleon scalar fields 
described in ^UJ; but disqualifies the construction as something that could change large-scale 
structure formation. (Note, for instance, that even a very low concentration of dark matter 
particles of a particular scalar charge would screen the scalar force to uninterestingly small 
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length scales via the chameleon mechanism [TT]). 

Before the D-term potential © spoils the story, what makes the model (J3J) work as a 
model that includes screening particles is that the space of supersymmetric vacua is not a 
manifold: three branches join at the point where <pi = <p2 = 4>s = 0. Many supersymmet- 
ric field theories exhibit a similar feature of multi-branched moduli spaces, and a number 
arise from string theory constructions. We will return to examples in this general class in 
section 12.41 

2.2 Supersymmetry breaking 

If supersymmetry is broken, then presumably a potential for otherwise massless scalars 
is generated at the soft breaking scale. Following conventional notions of supersymmetry 
breaking, one would predict that this scale could be as low as about 0.1 eV — some 28 orders 
of magnitude below the Planck scale. But for the scalar force to have a range on the order 
of a megaparsec, we need the scalar mass to be some 28 orders of magnitude smaller still! 
Thus, in considering supersymmetric theories with flat directions as candidate ingredients 
in a dark sector with scalar forces, we are entertaining the notion that supersymmetry is 
unbroken in the dark sector, or at least that the scalar somehow remains very light. This 
seems to be in clear conflict with standard notions of supersymmetry breaking, which hinge 
in part on the concept of field theoretic naturalness. 

On one hand, in light of the cosmological constant problem, we are disinclined to rule 
out a construction solely on the basis of difficulties that arise after supersymmetry is broken 
- particularly when it provides an interesting class of alternatives to ACDM. 

On the other hand, it is intriguing to note that in AdS / CFT [F2\ E3 E] , scalar fields 
whose masses are comparable to the scale of spacetime curvature are dual to operators in 
the field theory whose dimension is order unity. The precise relation for AdS^/CFT^ is 
m 2 L 2 = A(A — 4) ^2 E], where L is the radius of curvature, i.e. the analog of the Hubble 
length in the static AdS$ geometry. AdS/CFT is a general non-perturbative equivalence 
whose validity does not rely on supersymmetry. Were it not for this duality, the presence of 
scalar fields with mass ~ 1/L (even after quantum corrections) would be very surprising in 
a limit where L is much larger than the Planck scale. With the duality in hand, such scalar 
fields are precisely as natural as scalar operators whose dimension remains finite in the large 
N limit that corresponds to L ^ 1/Mp\. 

It has been suggested [T^] that AdS / CFT may be extended to a correspondence between 
dSi and a three-dimensional Euclidean quantum field theory. In such a correspondence, 
scalar operators in the field theory with finite dimension would correspond precisely to 
scalars whose masses are on the order of the Hubble scale — the main ingredient that we 
need for our story. 4 Eventually, this welcome conclusion would have to be reconciled with 
the loop calculations that support the field-theoretic notions of naturalness. For now, we 

4 This line of argument, starting with scalars in AdS, was suggested to us by C. Vafa. based on work with 
collaborators ^J. We thank him for his permission to review it here. 
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take it as a good reason to believe that light scalars are indeed possible in a quantum theory 
of gravity. 5 



2.3 The Kahler potential 

Even before supersymmetry breaking, it is possible to run into trouble with coupling to 
the visible sector. The difficulty is the Kahler potential K, which can receive quantum 
corrections at any order in perturbation theory. 

Let us recall why it is important to have some control over the Kahler potential 
Yukawa couplings in the component lagrangian are accompanied by a factor e K ^ 2 , so quark 
masses have some dependence on scalars in the dark sector. The proton mass will also 
depend on such scalars: if 

K = (constant) + Re(£0)/M PZ - \<j)\ 2 /M Pl + • • • , ( 7 ) 
where £ is a dimensionless complex number, then 



m v ~ m p + e p 



Re£0 |0| 2 , (Re£0) 2 



+ 



2Mpi M Pl 8M Pl 



where rh p is the proton mass when = 0, and e p is of order the u and d quark masses, 
that is, approximately 7MeV. In order not to violate current experimental constraints on 
violations of the Equivalence Principle, it is necessary to have both \<f>/Mpi\ ^2x 10 -5 and 
|£| ^ 2 x 10~ 5 . The first is not a problem based on estimates presented in the deviation of 
<f> from its background value is plausibly less than 10 -6 Mpi. The second appears to be more 
difficult to ensure: it amounts to the statement that just when extra massless modes appear, 
there is no linear term in the Kahler potential. If <fi is charged under a gauge symmetry, then 
this linear term is impossible — but then one has the D-term potential associated with the 
gauge coupling to contend with. It is perhaps illustrative to consider the case of the cone 
over SU(3) /U(l) 2 . There, at the classical level, the Kahler potential is K = — 3log V where 
V is the volume of the Gi manifold and <fi is one of the scalars <pi in a regime where the 
singularity is resolved enough to avoid Planckian curvatures. The cone volume is maximized 
when it is singular, leading us to expect K ~ constant — $^$: this is confirmed by the 
absence of the linear term because of gauge symmetries. 



2.4 Models with a non-abelian gauge interaction 

Let us now describe supersymmetric field theories generalizing (0) which avoid the pitfalls 
of D-terms and linear terms in the Kahler potential described in previous subsections. 

5 The loophole that permits light axions would imply rather special forms for the higher-point correlators 
of the dual operators. The above argument could be more precisely phrased as saying that light scalars with 
non-derivative couplings are dual to a class of operators whose dimensions should not all become large in a 
large N limit. 
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The combination of D-term and F-term potentials can leave protected flat directions in 
the space of supersymmetric vacua. As an example, consider SU (2) gauge theory with one 
adjoint flavor X a b and several fundamental and anti-fundamental flavors, qf and qi a . The 
renormalizable superpotential 6 is 

W = y s ^ a X a b q^ + ^Mitfiaq? + l -m x X\X b a . (9) 

If Mij = = mx, then the F-term potential vanishes either when X = and q^ a q\ = 
(referred to as the Higgs branch because q and q are allowed to be non-zero, which completely 
breaks the gauge symmetry) , or when q = = q (referred to as the Coulomb branch because 
non-zero X breaks the gauge group to U(l)). The D-term potential adds quartic self- 
interactions on both branches, but flat directions remain: in particular, Vjy = tr[X, X^} 2 on 
the Coulomb branch, which (after accounting for gauge freedom) leaves a single-complex- 
dimensional moduli space. 

The dynamics of D3-branes near D7-branes is described by a field theory similar to @ : 
see for example [TZ| for details. It has M = 2 supersymmetry and moduli space composed 
of a Higgs and Coulomb branch that meet at the point where gauge symmetry is restored. 

Provided mx = and has zero eigenvalues, theories like Q are good candidates 
for supersymmetric extensions of (JIJ. The scalar is replaced by the Coulomb branch, 
and the screening particles are both the qf and qi a and the non-abelian gauge bosons that 
become massless when X = 0. These particles interact with each other, but this does not 
substantially change the dynamical screening phenomenon which we review in section |3] 
Furthermore, gauge invariance prevents linear terms in X in the Kahler potential. And if 
has some eigenvalues near the Planck scale, the corresponding qf and qi a can play the 
role of massive dark matter. 

The above analysis is substantially unaffected by quantum loops provided the gauge 
interaction does not become strong in the infrared. This is guaranteed if there are enough 
light fundamental flavors $ and fy. 

The theory Q is evocative of M = 2 supersymmetric theories: with special choices of 
the parameters in (JUJ) — in particular, m x = — it becomes Af = 2 SU(2) gauge theory 
with several hypermultiplets in the fundamental representation. In the spirit of minimalism, 
it is interesting to inquire whether we can make do with the simplest such gauge theory, 
namely the case where there are no q and q fields. This theory has strong coupling in the 
infrared, but the dynamics is well understood due to (THJ. There are special points on the 
Coulomb branch where dyons become massless, but the SU(2) gauge symmetry is never 
restored (unlike the case when there are enough matter fields to make the gauge theory free 
in the infrared). The U(l) gauge boson is in a vector multiplet with a scalar a, where roughly 
X = aJs with J3 a generator of SU(2). At the point where magnetic monopoles become 

6 The superpotential @ is not the most general one allowed by renormalizability, but further generaliza- 
tions of it are not particularly interesting for our purposes. 
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massless, one may dualize to a magnetic description in which the light degrees of freedom are 
(in M = 1 language) Ajj, M, and M with a superpotential W oc MAjyM. It looks promising 
to regard as the scalar </> in and to let M and M be the screening particles. Other 
dyons might be the massive dark matter. But the Kahler potential — which is constrained 
by M = 2 supersymmetry — has a linear term in Ad, and this spoils the story as we have 
seen in section |2~^1 ' 

Finally, theories with M = 4 supersymmetry also have protected fiat directions, extra 
massless particles at enhanced symmetry points of spins 0, 1/2, and 1, and a protected Kahler 
potential. Moreover, they are realized on coincident D-brane world-volumes (in the absence 
of fluxes or other supersymmetry breaking effects). But explicit masses for the heavy dark 
matter particles would not respect the M = 4 supersymmetry. 

2.5 Summary 

Let us summarize: 

• If the dark sector is supersymmetric, massless scalars with a non-derivative coupling 
to non-relativistic dark matter become possible. 

• Supersymmetry breaking in the dark sector at the lowest natural scale, roughly 0.1 eV, 
tends to give scalars a mass on this order. Thus, naturalness arguments pose a problem 
for long-range scalar forces. 

• There are hints from AdS/CFT that scalars with mass comparable to the curvature 
scale are possible even without supersymmetry, despite the naturalness argument. 

• Absent an experimentally testable understanding of supersymmetry breaking, we pro- 
pose to search for viable models with supersymmetry, ultimately with a string theory 
origin. Viable means that the model shouldn't have problems before supersymmetry 
breaking. 

• One such problem is the lifting of flat directions by D-terms when a global symmetry 
becomes gauged: this is the fate of the W = y s <&i<&2<&3 model considered in section |2~T1 

• Another such problem is the presence of linear terms in the Kahler potential at the 
point where extra massless particles arise: this is the fate of the model based on the 
Seiberg-Witten solution considered at the end of subsection 12.41 

• Both problems are avoided by a class of models, considered in subsection 12.41 with a 
cubic superpotential and an adjoint chiral superfield which provides the flat directions. 

7 The offending term in the Kahler potential arises because there is a linear term in the prepotcntial: 
a = dJ-o/dao ~ do + —cld loga^ with nonzero oq, and the Kahler potential is IvadTo / dA£)A^ D . 
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• Models of this general type describe the dynamics of D3-branes near D7-branes, so 
there is good reason to expect that they arise in quasi- realistic compactifications of 
string theory. 

One may reasonably ask whether supersymmetry and/or string theory are necessary to the 
discussion. Indeed, one may take a completely phenomenological approach; yet it seems 
that string theory provides just the right ingredients for dynamically screened scalar forces 
acting with a strength comparable to gravity. Those ingredients are: moduli, supersymmetry, 
points of enhanced symmetry, and stable heavy objects such as wrapped branes. 

The bottom line is that string theory can help inform studies of an interesting class of 
generalizations of ACDM. The issue of supersymmetry breaking is not lightly to be dismissed. 
But we feel justified in bending the rules when the output is a class of models which may be 
tested against astrophysical data. 

3 Dynamical screening 

We turn now from the analyses in section |21 of the theoretical basis for the lagrangian (0) 
to a review and extension of the application of the lagrangian to the dynamical screening 
mechanism developed in pQ . In subsection 13.11 we recapitulate the basic mechanism. In 
subsection 13.21 we show that the screening particles can be excluded from halos formed from 
dark matter with a scalar charge. In subsection 13.31 we consider bounds on the various 
parameters of the model. 

3.1 The screening mechanism 

We assume that the propagating (source-free) part of may be neglected, and that the 
number density n s of the screening particles (the quanta of ^> s ) is large enough to drive 
close to zero almost everywhere, so the screening particles are relativistic. Since the source for 
— in the inhomogeneous distribution of the massive dark matter — varies with position on 
scales much larger than the Compton wavelength of the screening particles and much smaller 
than the Hubble length, the motion of a screening particle is close to adiabatic. That is, its 
energy e is conserved, apart from the effect of the general expansion of the universe. The 
screening particle mass is m s = y s \4>\ (assuming as usual that y s > 0), its velocity is v, and 
its conserved energy is e = m s / a/1 — v 2 . 

The equation of motion for 0, neglecting source terms from the non-relativistic species, 

is 

□0 = y 8 V s * s = y,n 8 (VT=v*) sgn = *^0 = -0 . (10) 

f T 

In the second equality of (jlOj) . one observes that ^ s ^ s = n s when v — 0, and that both ^ s ^> s 
and n s \/l — v 2 are scalars under Lorentz transformations. The next equality preserves the 
sign of 0, because the screening particles have positive energies. The characteristic mean 
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particle energy, from the average of the reciprocal Lorentz factor, is e s . The last expression 
defines a screening length, 

V Vs n s 

Now, operating in the quasi-static approximation — where we drop time derivatives of </> - 
let us re-introduce the two non-relativistic dark matter species, to get the wave equation 

V 2 = y+n+(t, r) - y_n_(t, r) + <J)/r 2 s , (12) 

where n± are number densities and y± are again positive. The first two terms on the right- 
hand side act as scalar charge densities and the last term produces the exponential cutoff in 
the scalar interaction of the dark matter. 

The spatial mean of the wave equation (|12|) gives the relation 

-y+n + + y_n_ 



^ ^ ±!1± ^ J ^ = ((vT3 ^ } ggn0) (13) 

y s n s 

where the inner brackets indicate averaging over the ensemble of screening particles to obtain 
a typical inverse Lorentz factor, while the outer brackets indicate averaging over space. We 
are assuming the magnitude of the ratio rj of scalar charge densities is less than unity, so 
\J\ — v 2 may be small almost everywhere. This means the screening particles are allowed 
to move so most remain relativistic, in positions where <fi is locked close to zero. Since the 
particles are conserved (apart from the small fraction captured in black holes) this condition 
is independent of redshift. One sees also that, as noted earlier, when the screenng particles 
are present the mean scalar charge density y+n+ — y~n~ need not vanish. 

It may be important that when <fi is close to zero nearly everywhere, so the screening 
particle nearly fill space, the general expansion of the universe causes e s to vary inversely 
as the scale factor a(t). This means the screening length r s increases linearly with a(t). 
Thus screening can occur at astronomically interesting length scales at low redshift but be 
insignificantly small at the epoch of formation of the observed anisotropy of the thermal 
cosmic microwave background radiation (the CBR). 

Let us consider finally the evolution of the peculiar velocities of the screening particles. 
When the particles are free to roam the general expansion of the universe stretches de Broglie 
wavelengths in proportion to the expansion parameter, a(t). Equation ()13J) says the particle 
mass is m s = \t]\e s , so the momentum of a particle varies as 



p = ye 2 - m 2 = eJl - (r]e s /e) 2 oc a(t) 1 . (14) 

If the bulk of the screening particles are relativistic to begin with, then their momenta and 
energies both scale with the expansion of the universe as a(t) _1 , the mass scales as a(t), and 
the screening particle peculiar velocities are independent of redshift. 
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3.2 Excluding screening particles from halos 



When clustering in the spatial distribution of the scalar charge density has grown large 
enough, (f) is pulled far enough from zero that y s \4>\ > e s . Where this happens the screening 
particles are excluded: their mass y s \<j>\ in this region is greater than their characteristic 
energy. Screening ceases in such a region. We find a condition for the screening particles 
to roam freely by using the solution to (|12|) valid in the limit where the length scale r v for 
variation of the source density y+n + — y-U- is much larger than the screening length r s . In 
this limit the left hand side of (|12|) may be neglected, and we have 

« r 2 s {-y+n + + j/_n_) . (15) 

Then the condition y s \(j)\ < e s with the definition (jlljlof r s is 

| | _ \y+n+(t,r)-y_n_(t,r)\ ^ ^ 

y s n s 

In the opposite limit, where the scale r v of clustering of the dark matter is small compared to 
r s , ()12|) becomes the Poisson equation, and the free roaming condition may be approximated 
as | r\ | £ (r s /r v ) 2 , that is, it allows a larger density contrast. 

At high redshift, where the dark matter distribution is close to homogeneous, (|16|) is 
satisfied everywhere if, as we are assuming, \rj\ in (|13|) is small. At low redshift, where 
the dark matter is strongly clustered, a region where n + is large will tend to drive out the 
\]/_ particles, so that y+n + ^> y_n_. Here we may write the condition that the screening 
particles roam freely in the presence of an overdensity of \& + particles with size r v 3> r s as 

^>^M = 1 + Mt,r), (17) 
y + n + n + 

where the number density n + (t,f) is expressed in terms of the mean density n+(t) and the 
contrast 5 + (t,r). If y s n s /y + n + = 2, for example, the screening particles may fill space in 
the early universe, when the mass density fluctuations are linear, and up to the time when 
the fluctuations start to become nonlinear on scales comparable to r s . If y s n s /y + n + = 100 
the screening would last into significantly nonlinear mass fluctuations, but would be broken 
today within the nominal virial radius r v ~ 300 kpc of a normal galaxy such as the Milky 
Way, unless r s 3> r v . 



3.3 Bounds on parameters 

Another condition to consider is that the relativistic energy density e s n s in the screening 
particles must not be large enough to spoil the standard model for the origin of prestellar 
helium and deuterium. We update the constraint in P under the assumption that the two 
non-relativistic dark matter components have the same mean densities of mass and scalar 
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charge. 

The result of expressing p s /r 2 s in terms of the Hubble parameter H and the energy 
fraction Q s in screening particles, with the assumption n + m + = n_m_ = p m /2, is 



fi s 35 ++ / y s n s 



(H r s ) 2 8 \y+n +i 



^ (18) 



where Q s and Q m are the density parameters in screening and the nonrelativistic ty± dark 
matter particles, and 

B ++ = si < 19 > 

is the ratio of the scalar and gravitational forces between two particles that are much 
closer than the screening length. The constraint from light element nucleosynthesis in the 
early universe allows relativistic mass beyond the standard model equivalent to about one 
two-component neutrino family [T^], or p s £ 0.23aT^, where the present CBR temperature 
is T Q = 2.725 K. At Hubble parameter H D = 70 km s _1 Mpc" 1 this condition is 

Q s < 1.1 x 10" 5 . (20) 

With dark matter density parameter Q m ~ 0.3 this bound in (|18|) is 

80 y+n+ , , 

^ < ^72 £ ^ M P C - (21) 
B+ 7 + y s n s 1 ' 

If B ++ = 1, so the scalar force is as strong as gravity at separations r < r S) and y+n + = 
0.01y s n s , which would keep the scalar pinned close to zero outside large galaxies, the screen- 
ing length could be r s ~ 1 Mpc, an interesting value for galaxy formation. 

If the screening particles were the decay products of massive nonrelativistic particles with 
half life much larger than one minute then the present mean energy density in screening 
particles could be larger than the energy density in the CBR without upsetting light element 
nucleosynthesis. But if the half life were shorter than 10 5 yr and p s were comparable to aT^ it 
would increase the cosmological expansion rate at decoupling, pushing the peak of the CBR 
anisotropy spectrum to smaller scales, and upsetting the present apparent concordance of 
constraints on the cosmological parameters. The constraint (J21J) thus seems to be reasonably 
general. 

To develop a more complete set of upper and lower bounds on the various parameters of 
the model, let us note that the equalities p s = n s e s and r 2 s = e s /y 2 n s can be combined to 
give (fTSj) and a second relation, 

3n s (H r s ) 2 = (-^r-) • (22) 
\ y s Mpi I 
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Let us start from a set of bounds intended to make the dark sector deviate in an acceptable 
but potentially interesting way from ACDM: 



V Tl 

> 1, 30 kpc < r s < 30 Mpc, - < B ++ < 10 . (23) 
y + n + 2 

The lower bound y s n s /y + n + > 1 ensures that if there is a scalar charge imbalance among the 
heavy particles over a very large region of the universe, the scalar field doesn't get pushed 
far away from zero. We could tolerate a much smaller value of this ratio if m + (4>) and 
m_(0) were perfectly linear, so that n + m + + n_m_ may be made independent of <fi: then 
an arbitrarily small n s locks <fi close to zero. But if departures from linearity in m±(0) arise 
from dimension 5 operators like (0 2 /Mpi)\i/ + \I/ + , then we do need y s n s /y + n + > 1 to lock 0. 

With the cosmological parameters in (f2T)|) and (|2"Tj) . equations (IT%|) and (J2"2"J) with (|2Tfy 
give 

V Tl 6 

■^-^ < 4000, tt s > 10~ 12 , 1.3 x 10" 11 < — s — < 9 x 10" 5 . (24) 

y+n+ y s M Pl 

A few comments are in order: 

• The bounds (|23|) are deliberately inclusive, and it may not be interesting to simulta- 
neously saturate the lower bound on r s and B ++ , as we have done, for example, to 
derive both bounds in (J2IJ). Thus in particular, y s n s /y + n + is constrained to be in a 
fairly narrow window. 

• As in j2j, we do not directly constrain y + or m+, but rather their ratio, entering the 
above bounds through the ratio B ++ of scalar and gravitational forces. 

• Likewise, we do not directly constrain y s or e s , but rather their ratio. 

If one assumes y s ~ 0(1), then the acceptable range (|2*H> for e s represents large energies: 
so large that at a redshift of 10 11 the typical particle's energy would be above the Planck 
mass. One way to arrange this is to suppose that the screening particles arise from decays of 
very heavy particles which have long but finite lifetime: for masses on the order of Mpi, the 
lifetime should be on the order of a month. An alternative — making e ~ 1 eV — requires 
n s ~ 0.1 cm -3 and y s ~ 10~ 20 , an exceedingly small value. 



4 The growth of structure 

We consider here the effect of the scalar interaction in two standard and simple models 
for structure formation: linear perturbation theory, which gives a good description of the 
growth of the clustering of mass and scalar charge at high redshift, and spherical symmetry, 
which illustrates effects of the nonlinear growth of structure. For the purpose of our con- 
siderations of the void phenomenon in the last subssection we will be particularly interested 
in the case where the two species of nonrelativistic dark matter have very different scalar 
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charge-to-mass ratios. The dominant theme in this case is that the light component quickly 
responds to any scalar field gradient, and the massive component more slowly responds to 
the resulting distribution of the light component. We discuss the possible relevance to the 
void phenomenon in subsection 14.31 This is not the only case of possible interest, of course. 
We briefly comment on other situations in sectional 

4.1 Linear perturbation theory 

The mean mass density in non-relativistic species is 

p m (t) = m + n + + m + n + + p b , (25) 

where the last term represents the baryons. Following standard practice we write n = n(l+S) 
for the spatially varying number density n of each species and its density contrast 5. The 
total density contrast S m in non-relativistic species satisfies 

Pm^m = rn + n + 5 + + m_n_5„ + p b 5 b . (26) 

In linear theory the condition \q\ < 1 in (|16|) means n s is close to homogeneous if it began 
that way, so we can take r s to be independent of position. We can drop the time derivatives 
of 4> because the inhomogeneities evolve slowly. Thus to determine </> we simply solve (112(1 . 

Following j2] , it is straightforward to check that in linear perturbation theory the Fourier 
modes of the number density fluctuations evolve according to 

5 q + 2H5 q = AnGp m ^ PpqfpSp ■ (27) 

p 

The dots represent derivatives with respect to proper time, so that Hubble's constant is 
H = a/a. The indices q and p run over the particles, the particles, and the baryons, 
and f p = n p m p / Y^ q n q m q is the mass fraction in species p. The dimensionless scale-dependent 
quantities (3 pq are 

where k is the wavenumber of the perturbation (so that k is dual to the coordinate x 
in the standard form of the cosmo logical line element, ds 2 = —dt 2 + a(t) 2 c?x 2 ), and the 
scalar charges are Q± = rky± for the dark matter and Qb = for the baryons. Equation 
(J2ZJ) neglects radiation drag and the pressure of the baryon gas, both of which are good 
approximations on the scale of galaxies after recombination and prior to galaxy formation. 
We have remarked that r s grows proportionally with a(t), that is, a/r s is constant, so (3 pq 
depends only on the species p, q and the comoving wavenumber k. The solutions are discussed 
in j2|: one introduces the matrix E pq = \J~Tpfipq\ffq-, which is real and symmetric, and whose 
eigenvalues £ c and eigenvectors c q determine perturbations A c = J2q c q\/ fq$q- When we can 
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neglect the energy density in radiation and the cosmological constant the modes vary with 
time as A c ~ t 2 W 3 , wne re 7 C± = (-1 ± ^1 + 24£ c )/4. 

This analysis readily generalizes to more than two species of dark matter, each with 
its own scalar charge-to-mass ratio, and to more than one scalar field, by adding to the 
dimensions of the matrices B pq (k) and E pq . The addition of gauge-mediated forces is easy, 
but, as discussed in |2j, it is not interesting because unlike gauge charges attract, as in 
an ordinary plasma, so charge separation is discouraged, just the opposite of the scalar 
interaction. 

It is instructive to specialize to the case where baryons are a negligible fraction of the 
total matter density, the cosmological constant A may be neglected, and there is overall 
charge neutrality in the non-relativistic species — that is, y+n + = y~n-, or equivalently 
f + y + /m + = f-y-/m-. Note that our assumptions imply /+ + /_ « 1. Then there is an 
adiabatic mode, where 5+ = <5_ = 5 m and £ = 1, which grows as 5 m ~ t 2 ' 3 (the usual CDM 
result); and there is an isocurvature mode, where 5 iso = 5 + — 5^ and 

* = 1 - P+- = " B +- = a r +V ~ w / 2 / 2 > ( 29 ) 

47rGm + m_ k z + cr/r^ 



which grows as 1+ V 1+24 £)/ 6 ; faster than the adiabatic mode if < 0, and slower if 
> 0. 

When radiation drag and matter pressure may be neglected the baryon density contrast 
obeys the equation 

5 b + 2H5 b = ATxGp m 5 m . (30) 

This standard result says that the baryons tend to approach the total mass distribution. 

If we do not impose the charge neutrality condition y+n + = y-n_ then adiabaticity is 
not preserved and the growing mode of the mass density contrast departs from the familiar 
t 2 / 3 behavior. To see this, consider the case where = and ignore the effects of r s and 
the baryons. For this choice of parameters the gravitational attraction and scalar repulsion 
of unlike particles just cancel: the two species of dark matter evolve with no effect on each 
other. In the growing mode, 5+ varies as t 2 " /+J ' 3 where 7 + = (—1 + y/\ + 24/3 ++ / + )/4, with a 
similar expression for 8-. Charge neutrality implies (3 ++ f + = 1, so in this case we recover the 
familiar t 2 / 3 growth law for the growing mode of the mass contrast, independent of the initial 
S±. Without charge neutrality, but still with j3 + _ = 0, an initially adiabatic perturbation 
evolves into a mass perturbation dominated by the contrast, 5+ or <5_, which has the larger 
7. If j3 + - < 0, then without charge neutrality an initially adiabatic perturbation may evolve 
into one in which 6+ and 5_ tend to have opposite signs. 

4.2 A spherical model for nonlinear clustering 

The spherical model offers a simple way to illustrate some of the effects of strongly nonlinear 
clustering. The numerical example presented in |2] shows the nonlinear growth of a mass 
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concentration out of an initially isocurvature (homogeneous mass density) perturbation with 
globally neutral scalar charge. Here we are interested in the nonlinear growth of charge 
separation. We focus on the case of two dark matter components, (+) and (— ), with very 

different charge-to-mass ratios, we assume B ++ is of order unity (so B ;§> 1), and we 

assume the length scales are small compared to r s so B is independent of scale. We present 
first some simple analytic considerations within the spherical model, and then a numerical 
example. 

Consider a neutral dark matter halo with initially identical distributions of the two dark 
matter components, in static equilibrium with no angular momentum. Since the scalar 
charge density vanishes there is no scalar force: the structure is the same as a standard cold 
dark matter halo. We will suppose almost all the halo mass M is in the (+) dark matter 
component. Now imagine varying the distribution of the light component while holding 
the massive component fixed. The variation of the energy of the light component under 
variations of its characteristic radius, R-, is dominated by the scalar field potential energy, 
U-, until the separation is large enough to make the kinetic energy of the light component 
important. When R_ is much less than the radius R + characteristic of the (+) matter the 
potential energy of the light component is negative, because it is dominated by the self- 
attraction of the (— ) matter. When R_ 3> R + the potential energy is positive, because it is 
dominated by the repulsion of the (+) matter that the (— ) matter sees as a central point-like 
charge. The limiting behavior thus is 

£l~- GM ' B++ at itL«i? + , U^~+ GM y ++ atitL>it! + . (31) 

One sees that there has to be an extremum of the potential energy at U- > 0. If, as we 
suspect, U- does not have a local minimum near the extremum, the (— ) matter is unstable 
against slipping off the halo of (+) matter. 

At central mass density p c the characteristic time for the light component to slip off the 
massive component is 

t- ~ {Gp c B ++ m + /m_)- l l 2 . (32) 

We are assuming m + ^> m_ and B ++ ~ 1, so t_ is much shorter than the dynamical time for 
the massive component, t + ~ [GpJ" 1 / 2 . This means a perturbation to the distribution of the 
light component grows by a considerable factor, ~ expt + /t_ ~ exp J m + / m_ ^> 1, before 
the heavy component can react. The conclusion is that, unless isocurvature fluctuations 
are very strongly suppressed, a well-mixed halo is a transient phenomenon: in spherical 
symmetry the light component will slide into a compact core, where it finds a balance with 
its kinetic energy, or else it will disperse to radii ^ r s , carrying with it the extra binding 
energy of the dark halo that ends up bound by the scalar force as well as gravity. 

We turn now to a numerical illustration of a variant of this effect, in the spherically 
symmetric development of charge separation during the growth of a dark halo out of a small 
initial density fluctuation. This solution assumes the halo is much smaller than r s , and it 
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ignores the mass in baryons. The halo is neutral, but there is a small initial charge separation, 
with contrast comparable to the initial mass density contrast. 

We simplify notation by scaling the particle masses and charges at fixed densities of mass 
and charge so y + = ?/_. This does not affect the dynamics, but it brings the expressions for 
the forces on a (+) particle caused by a (+) particle at distance r and by a (— ) particle at 
the same separation to the forms 

m + g ++ = (Gm 2 + + —j , m + g + ^ = -- \Gm + m„ - — j . (33) 

In a spherically symmetry, the inward accelerations of the (+) and (— ) dark matter compo- 
nents are 



( 777 

M±(r) 1 + ] + M T (r)(l - B) 



(34) 

m± J 

where the masses within physical radius r are M±(r). In a further simplification we have 
written B = B + _ = y 2 /(47rGm + m_) and, as before, we are assuming the system is small 
compared to r s . Equation (}34j) displays the two free parameters in the physics within our 
assumptions. To illustrate the situation under consideration we have chosen 

B = 6, m_/m+ = 0.01. (35) 



This corresponds to charge-to-mass ratios y j 'v4wGm + = 0.24 and y / \J A.irGm_ = 24. 

The remaining free parameters are in the initial conditions. The initial integrated mass 
contrast in the growing perturbation mode is 



_ M(< r) 3 2 fnr_ 

m<r ~ pV(<r) 5a x C ° S \2R 



(36) 



This refers to the mass within radius r; it is to be distinguished from the mass density 
contrast 5 m = 5p(r) / p. The function in the last expression vanishes at the outer radius R of 
the system, meaning the total mass is not perturbed from a background model that we take 
to be Einstein-de Sitter. The function is flat at r ~ R and r ~ 0, so the mass densities are 
close to homogeneous at the center and periphery. The larger density near the center requires 
that there is an intermediate zone where the mass density is smaller than the background 
model. This situation might approximate the formation of an isolated dark halo in a region 
with low mean density. 

The parameter a x in ()36|) is defined by the evolution of the mass density at r <^ R 
when the scalar force vanishes (or is canceled by identical distributions of the two species 
of charged particles). In this case, for a mass shell close to the center, a x is the ratio of the 
radius of the mass shell at the moment when the shell stops expanding to the radius at the 
initial conditions in equation ([36)1 . In our numerical example this central expansion factor 
is a T = 30. 
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Figure 3: A spherical model for the early development of a dark matter halo, comparing 
evolution in the standard cosmology — shown as a dashed line — and the interacting dark 
matter model with the parameters in equation (|35|). 

The initial integrated charge density contrasts in the two dark matter components are 



where the parameters f± = m±/(m + + m ) are the mass fractions in the two components. 
The initial mass contrast is S m<r = f + S + + consistent with (|36|). The initial charge 

contrast is 



The parameter representing a primeval charge imbalance is a = 0.5 in our example. This 
means the initial density contrast in the low mass negative component is 1.495 times the 
initial mass density contrast, and the initial density contrast in the positive species is 0.995 
times the initial mass contrast. 

In figure 01 the computed mass distributions are represented as the mean density p(r) of 
the mass M(< r) within a centered sphere of radius r, as in (J36j) . The dashed curve shows 
the standard model, with no scalar force, plotted at the moment when the mass near the 
center has stopped expanding and is about to collapse. Because mass shells are not crossing 
we have a simple analytic solution, which tells us that the center has stopped expanding at 
time t x = Znc?J 2 tij '4, when the central density is /Wier/Pouter = (37r/4) 2 times the density at 
edge, and where is the initial time in the Einstein-de Sitter model. 

The solid curve in the figure shows the effect of the scalar interaction on the mass dis- 
tribution. It is plotted at the same time t x as for the dashed curve, with the same initial 




(37) 
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conditions on the mass perturbation. We have put 1% of the mass in the component with 
large charge-to-mass ratio. This light matter responds to a scalar field gradient much more 
rapidly than does the heavy component. The larger initial central density in the light com- 
ponent thus quickly drives 1% of the mass into a tight central concentration. This central 
negative charge in turn pushes away the dominant mass component. In the figure there are 
no massive particles within the cusp in the solid curve at r ~ 0.45-R. Massive particles at 
the cusp are momentarily at rest. Further out, massive particles that initially were near 
the center are moving outward, and the rest are still falling in to form a dilute halo around 
the tight concentration of 1% of the mass. This parallels the behavior suggested by our 
discussion of fl3T}. 

If we had changed the sign of a in (jH7j) . to place an initial positive charge density near the 
center, the light particles would have rapidly moved away, leaving a positively charged dark 
halo that is more tightly bound than in the standard model because of the scalar addition 
to the gravitational attraction. Again, this parallels our discussion of (j31j) . 

We can conclude that, for the parameters under consideration here, very small primeval 
isocurvature fluctuations would ensure that the charges in a dark matter halo are always 
separated. Under spherical symmetry the separation is to a light compact core and dispersed 
mass or else to a near normal dark halo and a dispersed light component. A more realistic 
analysis would allow for separation in the manner of a squeezed watermellon seed, but the 
spherical model might be a useful approximation to the first generation of halos. If the 
initial charge imbalance in this generation were systematic, due to a small mean excess of 
charge in the light component, the spherical model would suggest there is some suppression 
of formation of the first generation of low mass dark halos. We comment on the possible 
observational significance in section 

4.3 The Void Phenomenon 

We consider here a possible relation between the situation we have been discussing, where 
most of the mass is in one of the dark matter components, and the observation that the nearby 
voids contain strikingly few galaxies of any kind. The observational situation is reviewed in 
It may be compared to the distribution of dark matter halos in the numerical N-body 
simulation shown in figure 1 in [20] • In the simulation there are relatively large dark matter 
halos that might be suitable homes for the normal galaxies — with masses comparable 
to the Milky Way — that contain most of the visible stars. These large halos appear in 
concentrations that resemble the observed clustering of normal galaxies. But the simulation 
shows that in the voids between the concentrations of large halos there are many relatively 
small halos that would seem to be suitable homes for dwarf and irregular galaxies. In the real 
world there are many more dwarf and irregular galaxies than normal ones, but the evidence 
is that they avoid the voids defined by the normal galaxies. This may not be the case in 
the ACDM cosmology: the simulations predict the presence of numerous potential homes for 
dwarfs in voids. Is there really a significant difference between the observed and predicted 
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relative distributions of dwarf and normal galaxies? If so, might the scalar force help explain 
the difference? 

The first question is put to the test in [20], who use the nearest neighbor statistic that is 
applied to the observations in PJ. The conclusion in [20] is that the predicted distribution 
shows little sensitivity to the halo mass, consistent with what is observed. The problem 
with this conclusion is that the typical distances in the simulation are an order of magnitude 
larger than in the observations, as one sees by comparing figure 6 in j2D] to figures 4 to 6 in 
[3]. Thus the meaning of the comparison is not clear. We conclude therefore that there is no 
established quantitative demonstration of an inconsistency between theory and observation, 
but there are distinct visual indications of a problem. 

If there is a problem there may be a resolution within the interacting dark matter model. 
Let us return to the linear perturbation theory discussed in section 5.1, and consider the 
case m + ^> m_ and B ++ ~ 1 . When the mass in the light component may be neglected, and 
we can also ignore the mass in baryons, the perturbation equations ()27|) and ()28|) become 



d 2 5 + n ad5 + . „_ r . . . _. 



d 2 5- hd5- . _ 



6+ + (<L - 6+) B ++ ^± 
m_ 



(39) 



The second equation says that the characteristic time for the evolution of <5_ is on the order 
of (GpB ++ m + /m„)~ 1 / 2 , as in (|32|). Since we are assuming m + ^> m_ and B ++ ~ 1 this 
is much shorter than the charactristic times for the growth of the mass contrast and the 
expansion of the universe, both of which are on the order of ~ {Gp)~ 1 ^ 2 . 

Under the conditions we are considering we are led to the following picture. On a time 
scale much shorter than the Hubble time the light component settles to minimize the sum 
of its potential and kinetic energy in the given distribution of the massive component. As 
the universe expands the kinetic energy in the light component is redshifted away, and 
this component tends to settle more deeply into its potential energy. The deepest minima of 
potential energy of the light component are in the regions of low mass density, the developing 
voids, where the scalar charge density of the massive component is lowest. The massive 
component responds to the clustering of the light component on a time scale comparable 
to the Hubble time. The concentration of the light component in the protovoids promotes 
evacuation of dark halos from the voids. In the ACDM cosmology dark halos tend to leave 
low density regions, as part of the general growth of clustering. The problem is that the 
tendency seems to be weaker than wanted. The scalar force would enhance this tendency, 
and so perhaps resolve the void problem. 

Baryons bound to dark halos in protovoids would be encouraged to leave with the massive 
dark matter. The only difference from ACDM in the behavior of baryons not bound to dark 
halos arises from the lower mass density in dark matter in voids predicted by our model - 
a small effect. It may encouraging therefore that low surface density HI clouds are observed 
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in nearby voids 

5 Discussion 

There is little doubt that a number of constructions exist in supersymmetric field theory 
and string theory which give rise to scalar forces screened by light particles, as in (JTJ. 
Incorporating non-relativistic dark matter particles coupled appropriate to the scalar seems 
unlikely to pose a difficulty, given that wrapped branes have roughly the right properties. 
Furthermore, it seems likely that the pitfall of linear terms in the Kahler potential can be 
avoided. What is clearly more difficult is to build such a construction into a compactification 
which also includes the Standard Model and which stabilizes those moduli which couple 
directly to the visible sector. Most difficult of all is to see in detail how scalar forces with 
sufficient range and appropriate couplings to affect large scale structure formation survive 
supersymmetry breaking. 

Notwithstanding these challenges, it appears that string theory provides at least the 
right ingredients for models that would include interesting scalar forces in the dark sector. 
Thus we propose a refinement of the moduli-fixing program in string theory: rather than 
requiring the elimination of all moduli, the objective should be to stabilize precisely those 
which couple to the visible sector in such a way that the absence of adequate stabilization 
would violate experimental constraints. Then one should ask whether other moduli could 
be left unfixed by vacuum effects — only to be stabilized by the presence of some number 
density of relativistic and non-relativistic particles, possibly leading to an interesting variant 
on the standard cold dark matter cosmological model. If future observations decisively favor 
models of this type over ACDM, then astronomy may become a window into string physics. 

We have not provided an account of initial conditions for the variant of the cold dark 
matter model under discussion here. More specifically, we have not described a mechanism 
which will produce the heavy dark matter particles or the relativistic screening particles in 
the correct abundances, or set the typical energy e s of the screening particles in the rather 
high range needed if y s is of order unity. The primordial perturbation spectrum is less of a 
worry, since what we need is a predominantly adiabatic spectrum plus a slight isocurvature 
component, which can easily arise from inflation. And the most immediate issue, we believe, 
is the possible observational consequences or indications of a scalar force in the dark sector. 
The discussion in section 0] assumes most of the mass is in one of the dark matter components. 
We have remarked on two possible applications, to the development of small-scale structure 
before substantial separation of the light and heavy dark matter components, and to the 
development of voids as the two components become well separated. 

The illustration in figure |3] of the possible effect of the scalar force on the formation 
of the first generations of dark matter halos assumes the presence of screening particles, 
which allows the postulate that the mean scalar charge density differs from zero, with an 
accompanying screening length, r s , on the scalar force (as discussed in section 3). The 
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Jeans length for the baryonic matter after decoupling, A ~ 10 kpc, is within the bounds 
discussed in section 3, so we are allowed to consider observationally interesting values of 
r s . Our illustration assumes the light component has the larger mean charge density. It 
suggests that this produces early formation of compact concentrations of the light component 
and early suppression of formation of massive dark matter halos. Our illustration also 
assumes spherical symmetry, however, so these indications certainly will have to be checked 
by numerical simulations before we may consider applying them to the apparent excess 
production of low mass dark halos in the ACDM cosmology (|22]> [22]), or to the developing 
constraints from strong lensing of quasars on low mass halos in the halos of massive galaxies 
and in the field (as discussed in [21] and references therein). 

It should be noted that the astrophysical community has not been very concerned about 
the issue of dark halos in voids, because the relation between dark halos and the baryons that 
illuminate them is difficult to predict. Many dwarf galaxies are observed in regions where the 
ambient density is close to the cosmic mean, however, and it is reasonable therefore to ask 
why so few extreme dwarf galaxies are observed in nearby voids. The proposal here is that 
the scalar force has pushed the massive dark matter out of the voids. Since this scenario 
assumes near charge neutrality and adiabatic initial conditions the scalar force has little 
effect on the linear evolution of the large-scale structure probed by the measurements of the 
anisotropy of the CBR. The effect of the scalar force on smaller scales, which is represented 
by the parameter B ++ , must not be too strong, because there are normal- looking galaxies 
in low density regions. Also, the scalar force must not upset the evidence that the Lya 
forest, which would react to the scalar force only through its gravitational interaction with 
the dark matter, is a good tracer of the dark matter power spectrum (as discussed in [23] 
and references therein) . It will be interesting to see numerical explorations of these effects, 
under the assumption that r s £ 10 Mpc. 

We have not discussed here a different parameter choice, in which the oppositely charged 
dark matter components have near equal particle masses. We remarked in [2] that if m_ ~ 
m + then a close pair of galaxies with dark matter halos of opposite charge could be only 
loosely bound, the gravitational and scalar forces nearly canceling. This with suppression of 
the exchange of dark matter particles between oppositely charged halos would reduce stellar 
dynamical drag, increasing the lifetimes of oppositely charged binary galaxies. Such an 
effect has some apparent observational support in the abundance of binary galaxies. In the 
absence of the screening particles there would be no preferred separation of binary galaxies; 
the distribution of separations would be set by the process of their production. With the 
introduction of the screening particles one can choose parameters so a pair of oppositely 
charged dark halos finds an equilibrium separation at r ~ r s . There is no evidence of this 
preferred separation in the galaxy two-point correlation function, but a closer look at the 
observations might be worthwhile. 
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